Matrix Element Randomness, Entanglement, and Quantum Chaos 
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We demonstrate the connection between an operator's matrix element distribution and entangling 
power via numerical simulations of random, pseudo-random, and quantum chaotic operators. Creat- 
ing operators with a random distribution of matrix elements is more difficult than creating operators 
that reproduce other statistical properties of random matrices. Thus, operators that fulfill many 
random matrix statistical properties may not generate states of high multi-partite entanglement. 
To quantify the randomness of various statistical distributions and, by extension, entangling power, 
we use properties of interpolating ensembles that transition between integrable and random matrix 
ensembles. 
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Random matrices have long been used as statistical 
models for a host of complex systems in areas of physics 
ranging from quantum dots to field theory Origi- 
nally introduced by Wigner to predict the energy level 
spectrum of heavy atomic nuclei 0, random matrices 
were later applied to characterize quantum chaotic sys- 
tems Lit La • More recently, random matrices have been 
found essential in various protocols of quantum compu- 
tation and communication including quantum state su- 
perdense codingH], remote state preparation ||, data 
hiding schemes |?J, and entangled state production jij. 

Given the similarities between quantum chaotic oper- 
ators and random matrices, and the aptitude of random 
matrices for producing entangled states, several authors 
have proposed en tanglem ent production as an indicator 
of quantum chaos |lCtllll While there has been some 
debate 13], substantial numerical evidence suggests that 
quantum systems with underlying chaotic classical dy- 
namics produce entanglement at a faster rate than sys- 
tems with underlying regular dynamics. 

In this work we explore the distribution of unitary op- 
erator matrix element amplitudes, a statistical criterion 
of random matrices [Til ]. We show that the matrix el- 
ements are key to entanglement production. Further- 
more, we demonstrate that, in creating an operator, it 
is more difficult to reproduce the statistical properties 
of the elements of random matrices than other random 
matrix properties. Evidence is provided from matrix 
ensembles that transition between integrable and ran- 
dom matrices, known as the interpolating ensembles, and 
pseudo-random matrices. For these classes of matrices 
nearest neighbor eigenvalue spacing and eigenvector el- 
ement distributions converge to that of random matri- 
ces more quickly than the matrix element distribution. 
The entanglement production of these operators is corre- 
spondingly slower in converging to that of random matrix 



production. We also use interpolating ensembles proper- 
ties to quantify the randomness of a given statistical dis- 
tribution and, by extension, entanglement production. 
Finally, we discuss the role of the matrix element dis- 
tribution in the entanglement production of quantum 
chaotic operators: why the entanglement produced in- 
creases with time, and how the entanglement production 
differs from non-chaotic operators. 

The circular ensembles of random unitary matrices 
were introduced by Dyson J^j as alternatives to the 
Gaussian ensembles of random Hermitian matrices 0, 
llfij . As with the Gaussian ensembles, the circular en- 
sembles play a central role in modeling complex quan- 
tum systems. Unlike the Gaussian ensembles, however, 
elements of the unitaries are not independent random 
variables 0] and are thus more difficult to generate. 
Nevertheless, matrices of the circular unitary ensemble 
(CUE), the space of all unitary matrices, can be gen- 
erated by taking the eigenvectors of a Hermitian matrix 
belonging to the Gaussian unitary ensemble (GUE), mul- 
tiplying each by a random phase, and using the resulting 
vectors as the matrix columns |l4j . 

The squared modulus or amplitude of the CUE ma- 
trix elements follow a distribution equal to that of GUE 
eigenvector element amplitudes. Let c l k denote the kth 
component of the Ith GUE eigenvector. The distribution 
of amplitudes rj = \c l k \ 2 in the limit of infinite Hilbcrt 
space dimension, N — > oo, after rescaling to a unit mean 
is Pgue{u) — e~ v j where y — Nrj 0. Since rj is un- 
changed when multiplied by a phase, the distribution, 
Pcue{ x )i °f the rescaled amplitude of CUE matrix ele- 
ments x, is equal to Pgue{v)- 

CUE matrices can also be generated based on the Hur- 
witz parameterization using elementary unitary transfor- 
mations, E^'^((j),ip,x) with non-zero elements 



*To whom correspondence should be addressed; Electronic address 

weinsteiOdave.nrl.navy.miH 

^Electronic address: hellberg@dave.nrl.navy.mil 



E ki j) = 1. k = l,...,N, k + i,j 

£<P = e ^cos0, E% j) = e*sM 

E$ j) = -e-^sin^, E^ 3) =e-^'cos(f> (1) 
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to construct N — 1 composite rotations 
Ex = EW-Wfou^xt) 

E 2 = J B (A '- 2 ^- 1) (012,^12,O)^ A '- 1 ^(^O2,Vo2,X2) 

E N -l = E (la) (0JV-2.AT-1, 1pN-2,N-l, 0) X 
S( 2 ' 3 )(0jV-3,JV-l,^JV-3,iV-l,O) X 

••• ^^'^cw-i^cuv-i^-i) (2) 

and, finally, E/ cc/B = e la EiE 2 . . . En-l Angles tp, x, 
and a are drawn uniformly from the intervals 

< iprs < 2tt, < Xs < 2vr, < a < 2vr, (3) 

and (fi rs = sin _1 (£ rs 1 ^ 2r+2 - ) ), with £ rs drawn uniformly 
from to 1 [HI. Note that the 2 x 2 block E&'?n with 
m,n = i,j and r = is a random SU(2) rotation with 
respect to the Haar measure. 

An advantage of the above method is that it allows 
for a one-parameter interpolation between diagonal ma- 
trices with uniform, independently distributed elements 
and CUE This is done by drawing the above pa- 

rameters from constricted intervals 

< VVs < 2tt5, < Xs < 2tt<5, < a < 2ir5, (4) 

■ 1 /(2r+2) 

and setting <fi rs = sin ) with £ rs now drawn 

from to 5. The whole is multiplied by a diagonal ma- 
trix of random phases drawn uniformly from to 27r. The 
parameter 5 ranges from to 1 and provides a smooth 
transition between the diagonal circular Poisson ensem- 
ble (CPE) and CUE for nearest neighbor eigenvalue and 
eigenvector statistics 

The interpolating ensembles matrix element amplitude 
distribution, however, does not transition smoothly with 
6. Rather, the matrix element randomness lags behind 
that of the other criteria such that even large 5 induces 
matrix element distributions noticeably different from 
Pcue(x) as in figure ^ 

A practical measure of multi-partite entanglement is 
the average bipartite entanglement between each qubit 
and the rest of the system 0, [2(J , 

9 " 

o = 2--E 2V Wi. ( 5 ) 

i=i 

where pj is the reduced density matrix of qubit j. An 
operators' lack of matrix element randomness causes 
the distribution of Q, for computational basis states 
evolved under one iteration of the operator, to diverge 
sharply from Pcue{Q), the distribution expected for 
states evolved by CUE matrices. For the interpolating 
ensembles a high value of 5, such as S — .98, yields 
eigenvalue and eigenvector distributions extremely close 
to CUE. P(Q), however, is still far from Pcue(Q) due 




FIG. 1: (Color online) Distributions of nearest neighbor 
eigenvalue spacings (upper-left) and eigenvector element am- 
plitudes (upper-right) for matrices of the interpolating ensem- 
bles with 5 = .1 (dashed), .5 (dotted) and .9 (chained). The 
matrix element plot (lower-left) also includes the distribution 
for S = .98 (light solid line). For this 8 the eigenvalue and 
eigenvector distributions are indistinguishable from random 
(solid line), but the matrix element distribution converges 
much more slowly. This is manifest in the entanglement gen- 
erated by operating with 100 8-qubit matrices on all compu- 
tational basis states, shown on the lower right for 8 = .9 (o), 
and .98 (diamonds). 

to lack of randomness in the matrix elements as seen in 
figure 

The 8 characterized distributions of the interpolating 
ensembles provide a quantifiable randomness measure for 
not fully random distributions. Such a measure provides 
objectivity in reporting the distance a given distribution 
is from the random one and allows the randomness of 
different statistical criteria to be compared. Not all non- 
random statistics will follow one of these intermediate 
distributions, but we find them useful nevertheless. 

The randomness of an operator's matrix element dis- 
tribution can be quantified by fitting it with an interpo- 
lating ensemble eigenvector distribution. Since an opera- 
tors' matrix element distribution determines the amount 
of entanglement generated, reporting i5 of the eigenvec- 
tor distribution best fit provides a measure of entangling 
power. An operator or class of operators whose matrix 
element distribution conforms to a higher <5-eigenvector 
distribution produces, on average, more entanglement. 
To generate a semi-random state with a given amount 
of entanglement, a quantum computer programmer can 
apply an operator with the desired matrix distribution S. 
We note that the interpolating ensemble matrix elements 
themselves are not well estimated by the 5-eigenvalue dis- 
tributions. This most likely stems from constrictions of 
the Euler angles, Eq. (J2J, which also leads to a non- 
smooth entanglement distribution. 

To further explore the relationship between matrix el- 
ement amplitude distributions, entanglement, and other 
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statistical measures of randomness, we turn to pseudo- 
random operators [2ll I22I Eij ] . Despite being efficiently 
implcmcntablc on a quantum computer, pseudo-random 
operators have been shown to reproduce statistical prop- 
erties of random operators. The pseudo-random oper- 
ator algorithm is to apply to iterations of the n qubit 
gate: random SU(2) rotation to each qubit, then evolve 
the system via all nearest neighbor couplings [21]]. A 
random SU(2) rotation with respect to the Haar mea- 
sure is described by equations Eqs. (QJ and JjJJ. The 
nearest neighbor coupling operator we use is U nnc = 

e l ( 7r / 4 )X/j=i cr 2® CT J ^ where <r 3 z is the z-direction Pauli 
spin operator. The random rotations are different for 
each qubit and each iteration, but the coupling constant 
is always 7r/4 to maximize entanglement. After the m 
iterations, a final set of random rotations is applied. 

As to increases, the statistical properties of the pseudo- 
random operators converge to those of random operators 
|2l) . However, as seen in figure |3 these properties do 
not converge at the same rate. For n = 8, m — 2 the 
eigenvalue spacing distribution is well approximated by 
the 6 = .8 interpolating ensemble eigenvalue spacing dis- 
tribution, the eigenvectors by the 8 = .92 eigenvector ele- 
ment distribution, and the matrix elements by the 8 = .7 
eigenvector element distribution. Applying 100 m = 2 
maps to all computational basis states gives an average 
entanglement (Q) of .7004 compared to .9883 for CUE 
matrices. For m — 4 the eigenvalues are approximated 
by 8 = .9, while the eigenvectors are indistinguishable 
from the random distribution. The matrix elements are 
approximated by a S of .78 and (Q) = .8416. Operators 
of higher m lead to eigenvalue and eigenvector distribu- 
tions that are practically random. The matrix elements, 
however, lag behind with 8 = .88 and (Q) = .9339 for 
m = 8 and 8 = .98 and (Q) = .9790 for m = 16. 

The above examples illustrate the difficulty in generat- 
ing operators with a random distribution of matrix ele- 
ments. In both cases increasing one parameter, S for the 
interpolating ensembles and m for the pseudo-random 
operators, causes a convergence to CUE statistics. The 
convergence is relatively quick for eigenvalue and eigen- 
vector distributions but much slower for the matrix ele- 
ment distribution. The correspondingly slow convergence 
of entanglement production to Pcue{Q) is evidence of 
the connection between the matrix element distribution 
and entanglement generation. 

Quantum chaotic operators are known for their abil- 
ity to produce entanglement. Numerical simulations of 
two coupled subsystems demonstrate the greater entan- 
glement generation of chaotic versus regular quantum dy- 
namics [uj 0, and analytical results have been ob- 
tained through various methods 0,13 E3- Here, we are 
interested in a quantum chaotic operator's ability to pro- 
duce entanglement on a quantum computer. Applying a 
chaotic operator once will not, in general, produce entan- 
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FIG. 2: (Color online) Distribution of nearest neighbor eigen- 
value spacings (top left), eigenvector elements (top right), ma- 
trix elements (bottom left), and Q (bottom right) for 8-qubit 
pseudo-random maps of m = 2 (plus), 4 (x), 8 (square), and 
16 (o). The eigenvalue and eigenvector distributions converge 
to that of random matrices (solid lines) more quickly than 
the matrix element distribution. The eigenvalue spacing for 
m = 2 and 4 are fitted by the 8 = .8 (chained line) and .94 
(dashed line) eigenvalue distributions. The m — 2 eigenvec- 
tor distribution is approximated by S = .92 (chained line). 
The distribution of matrix elements for the m = 2, 4, 8 and 
16 pseudo-random operators can be approximated by 8 — .7 
(chained line), .78 (dashed line), .88 (solid line), and .98 (dot- 
ted line) respectively. The low entanglement production is a 
direct outgrowth of the matrix element randomness lag. To 
approach Pcue{Q) requires m ~ 40 |2l|| . 

glement on par with random operators 0, l2fl | . After a 
number of iterations, however, the average entanglement 
of initial computational basis states (Q(t)), where t is 
the number of operator iterations, can approach that of 
random operators. We show how this is attributable to 
the matrix element distribution. 

First, we note that an increase in (Q(t)) with time 
cannot be due to eigenvector or eigenvalue statistics. The 
eigenvectors of an operator are unchanged as a function of 
t and the eigenvalues become uncorrelated. The random- 
ness of the matrix elements of chaotic operators, however, 
increases with t. To demonstrate this assertion we revisit 
entanglement production of the quantum baker's map 
and explore other quantized chaotic maps. 

Initial computational basis states evolved under the 
quantum baker's map attain Q values close to the ran- 
dom matrix average only at large t 0- Figure |3p shows 
that this is due to the map's matrix element distribu- 
tion. The baker's map matrix elements do not at all 
resemble Pcue(x), however, for t = 100 the distribution 
is well approximated by the 5 = .9 interpolating ensem- 
ble eigenvalue distribution. For the n = 8 baker's map 
(0(1)) is only .3080, while (Q(100)) is .9597. More itera- 
tions lead to increased matrix element randomness which 
causes greater entanglement generation. 
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Other examples of quantized chaotic maps are 
the quantum sawtooth map [2?l Ejf . U saw = 
er^^rn? in ^(n- m f in ^ and thc quantum Harper 

map 0, U H = e iN '< cos ^' i / N ^e iN i cos ^P' N \ All ele- 
ments of the chaotic, k — 1.5, and regular, k = —1.5, 
sawtooth maps have equal amplitude. Operating on any 
computational basis state yields a state of Q — 1. For 
the chaotic sawtooth the matrix element randomness in- 
creases with t, such that at t = 50 the matrix element dis- 
tribution is practically Pcue(x) and (Q(50)) = .98826. 
For the regular sawtooth the entanglement oscillates 
wildly as seen in figure This stems from the lack of an 
asymptotic randomness for the matrix elements. 

The matrix elements for the chaotic Harper, 7=1, 
deviate only slightly from Pcue{x), and, thus, (Q(l)) = 
.9814. For t — 50 there is the expected increase in matrix 
element randomness and (Q(50)) = .9882. The regular 
Harper map, 7 = .1, matrix element distribution and av- 
erage entanglement also approach asymptotic limits as 
t increases. However, these limits fall short of the ran- 
dom matrix statistics. The matrix element distribution 
is well fit by the S = .9 interpolating ensemble eigenvec- 
tor element distribution and the average entanglement 
(Q(t — > 00)) ~ .95. We note that the matrix elements 
of the 100 times iterated baker's map and the m = 8 
pseudo-random operators also follow the S = .9 distribu- 
tion. Thus, (Q(50)) of the regular harper is very close to 
(Q) of these operators. 
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FIG. 3: (Color online) Average entanglement, Q, over all 8- 
qubit initial computational basis states as a function of time 
for quantum sawtooth maps, k = 1.5 (chained line) and k = 
— 1.5 (dotted line) and Harper maps, 7 = 1 (solid line) and 
7 = .1 (dashed line), compared to the random matrix average 
(horizontal dashed line). The chaotic maps quickly approach 
the random matrix average while the regular maps do not. 
The insets show matrix element distributions for the regular 
(light) and chaotic (dark) sawtooth maps at t = 50 (A), the 
regular (light) and chaotic (dark) Harper maps at t — 50 (B), 
and t = 1 (light) and 100 (dark) of the baker's map (C). The 
matrix elements for the t = 50 regular Harper map is well 
approximated by the S — .94 eigenvalue distribution. 



In conclusion, we have explored connections between a 
system's matrix element distribution and entanglement 
production. The interpolating ensembles demonstrate 
the difficulty in creating operators with randomly dis- 
tributed matrix elements and their statistical properties 
provide a useful measure of randomness. This measure is 
used to compare and contrast the randomness of statis- 
tical distributions from pseudo-random and chaotic op- 
erators. The connection to matrix element randomness 
provides a random matrix basis for increased entangle- 
ment production of chaotic systems as a function of time 
and the greater entanglement production for chaotic ver- 
sus non-chaotic systems. Though this need not always 
be the case, operators with similar matrix element dis- 
tributions appear, on average, to produce states with a 
similar amount of entanglement. 
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